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Abstract 

A previously published algorithm for trajectory tracking control of 
tethered wings, i.e. kites, is updated in light of recent experimental evi- 
dence. The algorithm is, furthermore, analyzed in the framework of delay 
differential equations. It is shown how the presence of system delay influ- 
ences the stability of the control system, and a methodology is derived for 
gain selection using the Lambert W function. The validity of the method- 
ology is demonstrated with simulation results. The analysis sheds light 
on previously poorly understood stability problems. 

1 Introduction 

During the past two decades researchers have taken an increasing interest in the 
industrial applications of tethered wings, i.e., kites. Two well-known applica- 
tions include the towing of ships [14] and wind- to electrical power conversion 
[13, 3]. An integral part of both technologies is autonomous kite flight along a 
prescribed trajectory, relative to the earth surface attachment point. Various 
control principles have been proposed, see, e.g, [8, 5, 2, 7, 16, 11]. This pa- 
per discusses the influence of system delays on trajectory tracking performance. 
The issue of system delays has been raised before in [8], where a wireless link 
is part of the control loop, and in [5], where actuator delay is identified. Up to 
now, however, no rigorous analysis of the effect of delay on system performance 
has been carried out. 

A small kite such as the one shown in Figure 1 is highly sensitive to control 
inputs, and hence to system delay. The ability to fly such small wings smoothly 
with minimal control expenditure therefore hinges on an accurate treatment 
of their fast dynamics. We use the theory of delay differential equations to 
incorporate system delay into the stability analysis of a kite trajectory tracking 
control algorithm. 

We base our analysis on the controller previously laid out in [2]. Since its 
initial publication, the algorithm has been applied in practice (Figure 1), and in 
the process, has undergone several revisions. In the next section, we commence 
by summarizing our algorithm in its most recent form. We then discuss the 
stability of an elementary linear differential equation with delay. Noting that 

*Baayen & Heinz GmbH, Sekr. ER2-1, HardenbergstraBe 36a, 10623 Berlin, Germany, 
jorn.baaycn@baayen-heinz.com. 



1 




our inner loop controller shapes the system into a delayed linear system of the 
same type, we proceed to identify its region of stability and, in a certain sense, 
optimal control gain. In the final section we demonstrate the applicability of 
the result in simulation. 



2 Trajectory tracking control 

In [2] we showed how the turning, or path, angle of the kite trajectory can 
be used to formulate a cascaded single-input single-output control problem. In 
this section we briefy recapitulate this result, and present several modifications 
which have turned out to improve performance in practice. 

First of all, the presence of the kite's tether constrains its movement to a 
sphere centered on the tether's earth surface attachment point. In our analysis, 
the radius of this sphere turns out to be irrelevant, and therefore we restrict our 
attention to the unit sphere centered on the origin. 

Restricted to the unit sphere, we can speak of geodesic, or great-circle, dis- 
tances 

dist(p, q) := acos(p • q), 

as well as of tangent planes. After defining a basis in the tangent bundle, we 
can measure the angle between any tangent vector and the primary basis vector. 
For vectors tangent to a curve, this angle is known as the turning angle [10] - 
in this paper denoted as 9 (Figure 2). 

In [2] we defined the steering control input, u, as the length difference of 
the kite steering lines. We showed how this control input u can be used to 
control the direction of flight, represented by the turning angle 9. Secondly, we 
used spherical geometry to derive an outer loop control law that selects a target 
turning angle, 9 t , such as to guarantee convergence to the target trajectory. 

A commonly prescribed target trajectory is the crosswind, lying figure of 
eight (Figure 1). This trajectory has a bounded winding number, thereby avoid- 
ing problems with lines winding around themselves and each other. 

We now proceed to summarize the inner and outer loop control laws. 
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Figure 2: The notion of turning, or path, angle 0. 



2.1 Inner loop 

In [5] it is, based on experimental evidence, claimed that 

6 = KV a u, (1) 

with proportionality factor K, airspeed V a , and steering control input u. 1 The 
linear structure of Equation 1 obsolctcs the incremental feedback-linearizing ap- 
proach previously taken by us in [2] , and thereby simplifies matters significantly. 

However, the airspeed V a is also a function of the turning angle, and hence 
also a function of the steering control input u. Implicit in Equation 1 is the 
assumption that the influence of the control input u acts on the turning rate 
9 more quickly than it acts on the airspeed V a . We return to this time-scale 
separation in Section 4. 

Proposition 2.1. Consider the system given by Equation 1. The control law 

ui=j^r(e t + p(0 t -e)), 

with gain P e R + , renders the Lyapunov function 

V:=\{6 t -ef 

decreasing along the flow. 

Proof. Expand dV/dt. □ 

Future work on the inner loop may focus on the incorporation of actuator 
dynamics using, e.g., backstepping [12]. 

1 In [5], the turning angle 8 is called the "orientation angle", denoted ip. 
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2.2 Outer loop 



As the kite moves along its trajectory, 7, the desired target point 7 t (s) must 
slide along on the target trajectory, j t . To this end, we continuously locate the 
target point closest to the present kite location, i.e., 

s(t) := argmin^ dist(7(i), 7t(«))- (2) 

This optimization problem is not convex. But, since we desire to track a 
smoothly changing target location, we welcome local minima. A simple gra- 
dient descent [9] performs well for this purpose. When initialized with the 
solution from the previous timestep, the method converges within a handful of 
iterations. Unlike the method we originally proposed in [2], this method does 
not suffer from drift of the target location. 

In the following we will, for brevity, identify 7 t with the composition 7 t o s. 
Once a target point has been determined using Equation 2, a target turning 
angle needs to be selected. To this end, let T t denote the unit vector tangent 
to the target trajectory. Let T :— R T T t , where the operator R rotates 7 to 74 
along their connecting geodesic, so that T is an element of the tangent plane 
at 7. Let t := 7 x T, so that the vectors T and t form a basis of the tangent 
plane at 7. Finally, denote 

ii ■= it - {it ■ 7)7- 

By selecting the target flight direction as an appropriately weighted sum of 
the "direction along the target" T, and the "direction towards the target" t, we 
obtain tracking convergence: 

Proposition 2.2 (Baayen, 2012 [2]). The turning angle 9 t determined by 

cos6* t ei + sin# t e 2 = normalizc(T + L(t ■ 7^)t), (3) 
with gain L £ R+ , renders the geodesic distance 

W := dist(7,7 t ) 

decreasing along the flow. 2 



3 A delay differential equation 

We now introduce an elementary linear delay differential equation, which in 
Section 4 will prove to be a valuable tool for analyzing the stability of our kite 
control system in face of delay. 
Let 

§(t) = *y(t-r), (4) 

with k £ R, a linear differential equation with delay r. With the change of 
variables s = t/r and x(s) = y(rs), we obtain the equation in its canonical form 

^(s) = ax(s - 1), (5) 



2 In the original paper, the required normalization is missing. For practical applications 
this does not matter, as the trigonometric computation of 8t is only a function of the direction 
of the vector, not of its norm. Thanks to Dipl.-Ing. Claudius Jehle for pointing this out. 
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with a — rk. 

Unlike its relative without delay, Equation 5 admits oscillatory solutions. 
For example, when a = —n/2, it is readily verified that 




provides such a solution. Hence it is not the case that a negative coefficient in 
Equations 4 or 5 guarantees convergence to the zero-solution. We must therefore 
consider these equations in more detail. 

It turns out that the solution space of Equation 5 is infinite-dimensional 
[6]. This circumstance is hardly surprising, considering that we need to specify 
initial conditions over an interval. Any solution can be written as 

n 

with the coefficients A„ determined by the initial conditions, and a n the roots 
of the characteristic equation 

a - aer a = 0. 

The roots of the characteristic equation, a n , are given by the set-valued Lambert 
W function: 

K} = W(a). 

The real part of the first several branches of the Lambert W function is shown 
in Figure 3. 




a 



Figure 3: Real part of the Lambert W function, branches —5 to 5, with the 
principal branch highlighted in bold. 

We proceed to analyze the stability of Equation 5. Let Wk denote the kth 
branch of the Lambert W function. Figure 3 hints at the following Lemma: 
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Lemma 3.1 (Shinozaki and Mori, 2006 [15]). The real part of the principal 
branch of the Lambert W function is no less than that of any other branch, i.e., 

max{3fi(VK fe (z))|/c = 0, ±1, ±00} = 5ft(W (z)). 

Lemma 3.2. The zero-solution of Equation 5 is globally stable if 

-tt/2 < a < 0. 

Proof. The proof follows directly from Lemma 3.1, Figure 3, and the observa- 
tions that W(-7r/2) 3 itt/2 and W(0) 3 0. □ 

Lemma 3.3. The quickest worst-case convergence to the zero-solution of Equa- 
tion 5 is guaranteed if 

1 

a = , 

e 

irrespective of the initial conditions. 

Proof. The proof follows directly from Lemma 3.1, Figure 3, and from the fact 
that W is not differentiable at — 1/e [4]. □ 



4 Trajectory tracking control with delay 

In the present section we modify the results from Section 2 to incorporate a 
system delay r. 

First of all, we recall the time-scale separation of the turning rate and the 
airspeed V a implicit in Equation 1. The turning angle 9 specifies the direction 
of flight. Considering that the direction of flight acts on the position indirectly 
through integration, we can expect another time-scale separation between the 
flight direction and the position dynamics. 

In the following we show how this structure can be exploited to analyze 
stability in face of system delays. 



4.1 Inner loop 

The outer loop control law given by Equation 3 is a function of position only. 
By time-scale separation, this implies that the resulting target turning angle 9 t 
varies slowly compared to the time constant of the inner loop. In the present 
discussion of the inner loop we therefore assume that dO t /dt = 0. 

Exploiting this assumption and following Section 3, we obtain the following 
stability criterion: 

Proposition 4.1. Consider the turning rate law with delay t: 

e(t)=KV a (t)u(t-T), (6) 

cf. Equation 1. Assume that d0t/dt — 0, and let P e (0, ^). Then the control 
law 

renders 9 — > 9 t . Furthermore, this convergence is optimal in the sense of Lemma 
3.3 when 
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Proof. Substituting Equation 7 into the system given by Equation 6, and using 
our assumption that d6 t /dt — 0, we obtain the controlled system 

±(0 t -e)(t) = -P(6 t -6)(t-T). (8) 

From Lemma 3.2, we obtain that the controlled system given by Equation 8 is 
stable if -tt/2 < -Pt < 0. □ 

This quantifies what we already know from practice, namely, that the inner 
loop control gain P cannot be arbitrarily large. Proposition 4.1 aids us in 
selecting a gain appropriate for the specific system implementation at hand. 

4.2 Outer loop 

The outer loop treats the target direction 9 t as a virtual control input, capable 
of being commanded pseudo-instantancously. This assumption is motivated by 
the time-scale separation. The target turning angle 9 t will be commanded to 
the inner loop without delay within the same control software cycle. 

Furthermore, we assume the system delay r to be insignificant compared to 
the time scale of the position dynamics. As we will be able to verify with the 
simulation results from the next section, the system delay is typically a small 
fraction of the flight time of a figure of eight. 3 Hence, we do not include delay 
in our analysis of the outer loop. 

5 Simulation results 

In [1] we show how one may simulate a kite system using potential flow com- 
putations and a discretized tether. Using an identical system setup, we now 
present several simulation results. 

We start out by presenting a reference result without delay. Selecting an 
outer loop gain of L — 5 and an inner loop gain of P = 10, the target trajectory 
is tracked, for all practical purposes, perfectly. See Figures 4 and 5. 

By introducing a delay of r — 0.2 s, the gain P = 10 lies outside of the range 
given by Proposition 4.1. As expected, this is reflected in the instability of the 
turning angle. See Figures 6 and 7. 

Finally, by choosing L = 5 and the optimal inner loop gain P = l/(er), we 
obtain the convergent trajectory and stable turning angle tracking of Figures 8 
and 9. We note, however, that trajectory tracking performance suffers somewhat 
at the expense of stable turning angle tracking. 

We conclude that Proposition 4.1 provides a useful tool for understanding, 
and eliminating, delay-induced instabilities observed in kite control systems. 

6 Concluding remarks 

In this paper we have presented, based on recent experimental evidence, a sim- 
plification of our previously published cascaded kite trajectory tracking control 
algorithm. Furthermore, we have shown how the Lambert W function can be 

3 In contrast, a kite can be made to spin rapidly about its tether. 
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Figure 4: Trajectory with t = s, Figure 5: Turning angle (solid) with 
L = 5, and P = 10. target (dashed), r = s, L = 5, and 

P = 10. 




Figure 6: Trajectory with r = 0.2 s, 
L = 5, and P= 10. 




Figure 8: Trajectory with r = 0.2 s, 
L = 5, and P = l/(er). 
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Figure 7: Turning angle (solid) with 
target (dashed), r = 0.2 s, L = 5, and 
P = 10. 
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Figure 9: Turning angle (solid) with 
target (dashed), r = 0.2 s, L = 5, and 
P = l/(er). 
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used to analyze the stability of kite trajectory tracking control systems with 
delay. We demonstrated how such an analysis leads to a methodology for gain 
selection in face of delay. Finally, we demonstrated the validity of our approach 
with simulation results. 

Future research may focus on the influence of actuator dynamics on inner 
control loop performance, and similarly, on the influence of inner loop dynamics 
on the performance of the outer loop. 
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